Acceptor and donor impurities in GaN nanocrystals 
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We investigate acceptor and donor states in GaN nanocrystals doped with a single substitutional 
impurity. Quantum dots (QD's) of zinc- blende structure and spherical shape are considered with 
the radius ranging from 4.5 to 67.7 A. The size-dependent energy spectra are calculated within 
the sp 3 d 5 s* tight-binding model, which yields a good agreement with the confinement-induced blue 
shifts observed in undoped QD's. The computed binding energy is strongly enhanced with respect 
to the experimental bulk value when the dopant is placed at the center of the smallest QD's. It 
decreases with increasing QD size following a scaling law that extrapolates to the bulk limit. In 
order to estimate the degree of localization of the bound carriers we analyze their wave functions 
and average radii. The resulting physical picture points to a highly localized acceptor hole, mostly 
distributed over the nearest-neighbor anion shell, and a much more extended donor electron. We also 
study off-center impurities in intermediate-size QD's. The acceptor binding energy is approximately 
independent of the dopant position unless it is placed within a surface shell of thickness of the order 
of the bulk Bohr radius, where the ionization energy abruptly drops. On the contrary, the donor 
binding energy gradually decreases as the impurity is moved away from the center toward the QD 
surface. 



I. INTRODUCTION 

Group-Ill nitride semiconductors such as A1N, GaN 
and InN have drawn considerable attention in recent 
years because they generate efficient luminescence and 
their band gaps include the entire visible spectrum and 
a part of the ultraviolet (UV) . As a result of intense re- 
search on these compounds, visible light-emitting diodes 
and UV laser diodes are available nowadays, opening the 
way to further developments in optoelectronics, specially 
in optical storage technology. These materials are also 
advantageous in high-power and high-temperature elec- 
tronics. We focus on GaN, since it is the most outstand- 
ing nitride of this group, widely used and researched^— 

The great interest in semiconductor quantum dots 
(QD's) is based on the possibility of tailoring their elec- 
trooptical properties through the control of size and dop- 
ing. Nanocrystalline~~ 7 and self-assembled^— QD's of 
GaN have been synthesized and show the confinement- 
induced blue shift. Recently, Eu-doped GaN nanocrys- 
tals (NC's) yielding a bright and narrow Eu 3+ ion flu- 
orescence line have been shown to be useful as an effi- 
cient marker of important biological processes^ 4 How- 
ever, to our knowledge, no experimental study of donor- 
or acceptor-doped GaN QD's has yet been reported. 
Theoretically, spherical NC's containing a substitutional 
impurity have been investigated previously* 1 ^— Both 
donors and acceptors have been studied within a vari- 
ational approach (VA)2£ and the effective-mass approx- 
imation (EMA) jis - 1 V 1 - 27 , 29 Acceptors have also been 
treated in the atomistic framework: Effective bond- 
orbital^ and tight-binding (TB) 2? models. 

The electronic structure of a doped QD is highly sen- 
sitive to the position of the impurity atom. But most of 
the theoretical studies reported so far have focused on 
on-center impurities. The impurity position dependence 



has been investigated only for donors in the continuous 
models, VA and EMA, which are known to be inadequate 
for the electronic structure of small-size QD's and highly 
localized acceptor states even in the bulk. In the present 
work we report a TB model calculation of the electronic 
structure for both on- and off-center donor and acceptor 
impurities in GaN NC's of diameter up to 13.5 nm. A 
preliminary study reported previously 2 ^ was limited to 
only on-center acceptor impurities in small-size NC's up 
to 4.5 nm in diameter. Also, it was based on the simple 
semi-empirical sp 3 s* model of Vogl et al£Q In contrast, 
the calculations presented here are based on the extended 
sp 3 d 5 s* model of Jancu et al.^ which provides a more 
satisfactory description of the band structure in the bulk. 
For the acceptor impurity we focus on Mg as an exam- 
ple. It substitutes for a Ga cation in the zinc-blende lat- 
tice. The hole binding energy in the bulk is known to be 
236 meVi^ 2 . For the donor we consider O, which replaces a 
N anion yielding an electron binding energy of 33 meV in 
the bulki^ The paper is organized as follows: In Section 
II we present an outline of the theoretical model; in Sec- 
tions III and IV we show the binding energy and average 
radius of the bound carrier as a function of the dopant 
location and the QD size; in Section V we conclude by 
summarizing the main results. 



II. THEORY 



The GaN host material is described by a semi-empirical 
sp 3 d 5 s* TB Hamiltonian which includes the spin-orbit 



interaction, and was first introduced by Jancu et 
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where i and j denote the atomic sites, /i and v the single- 
electron orbitals, and h % ^ v the hopping and on-site ener- 
gies; Xi is a constant standing for the spin-orbit interac- 
tion. The tight-binding parameters account for the elec- 
tronic structure of bulk GaNi22 

We need to take into consideration the dangling bonds 
on the surface of the quasi-spherical crystallite, since they 
introduce states within the band gap. The simplest way 
to move these states several eV's away from the band gap 
region consists in binding every dangling bond to a hy- 
drogen atom^ the Ga-H and N-H distances are taken as 
one half the respective covalent bond lengths. This sim- 
ple passivation procedure is more detailed in a previous 
work'' 4 and references therein. 

The effective Bohr radius a B of the hole bound to an 
acceptor in bulk GaN is typically small: In the case 
of Mg (binding energy Eg = 236 meV) it is 4.5 A. 
The Bohr radius a B of the electron bound to a donor 
is much larger: a d B = 27.7 A for O (Ef = 33 meV). 
These values of the Bohr radii arc obta ined in the hydro- 
genic model: ob = S/ y / 2Ef ) (bulk) m*, where the effec- 
tive mass is 0.8 mo for the bound hole? 8 - and 0.15 mo for 
the bound electron^ Since the donor Bohr radius is over 
six times larger than that of the acceptor, we need to 
construct large enough QD's to follow the confinement 
effects on donors. The resulting Hamiltonian matrices 
are huge and sparse. For example, our largest QD con- 
sists of ~ 113, 000 atoms and gives a matrix of size over 
2 • 10 6 x 2 • 10 6 . Furthermore, the new (lower) symme- 
tries introduced by the off-center impurities are difficult 
to handle within a group-theoretical framework, as in 
Ref. H|| Hence, we resort to Lanczos algorithms^— to 
calculate the few dopant levels we are interested in. The 
numerical libraries that we use are freely available at the 
Arpack site^ 

The impurity is modelled by a screened Coulomb po- 
tential: 



given by the simplified Penn model 4 ^: 
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where is the distance between the lattice site i and the 
impurity, R the NC radius, and e(r i; i?) the dielectric 
function within the crystal. Here Uq stands for the cen- 
tral cell potential which obviously cannot be described 
in terms of the long-range Coulomb potential that rep- 
resents the impurity atom in its crystalline environment 
as an effective point charge. In the absence of any pre- 
cise theoretical model, as usual, we shall treat Uq as 
a phenomenological parameter. Position-dependent di- 
electric functions have been proposed by Resta 3 -^ and 
by Franceschetti and Troparevsky.— The position depen- 
dence is, however, restricted to a small volume around the 
dopant, typically smaller than the distance between two 
nearest-neighbor atoms, 3 ^ beyond which the permitivity 
recovers the long-distance behavior. When restricted to 
QD's that means simply a value dependent on the QD 
size. We then take a size-dependent dielectric function 



e(R) = 1 



sb - 1 
1 + (£/Rp ■ 
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The parameters entering this model are 7 = 2, the 
material-dependent £ = txEf /kpE g = 7.43, by using 
E g = 3.3 eV and E F = 16 eV (Ref. Q), and the bulk 
dielectric constant^ £s(0) = 8.9. 

The acceptor and donor binding energies are respec- 
tively defined a a 28 ' 43 



E n — El 



HOMO, 



Ef; = E WM o - E d , 



(4) 



where -Ehomo (R) is the energy of the valence band-edge 
state, the so-called highest occupied molecular orbital 
(HOMO), and -Elumo(-R) the energy of the conduction 
band-edge state, the lowest unoccupied molecular orbital 
(LUMO), in the undoped NC. E a (R) and E d (R) corre- 
spond to the same states, but in the doped NC. 

We also compute the average radius of the bound car- 
rier, r(i?), which is related to the effective Bohr radius 
at large R: 

-/n \ 3a,B 1 

r(R^ oo) = cc - . (5) 

2 ^ (bulk) 



III. RESULTS: IMPURITY AT THE 
NANOCRYSTAL CENTER 

We first study the electronic structure of undoped 
QD's. The calculated energy levels Ehomo(R) and 
Elxjmo(R) are plotted against the NC radius in Fig. Q] 
They are the reference levels for calculating the impu- 
rity binding energies (see Eq. (Q}). Note the dramatic 
increase of the band gap with decreasing QD size. This 
confinement-induced blue shift is in accord with the re- 
ported spectroscopic studies of GaN NC's.— Quanti- 
tatively, we obtain a satisfactory agreement with exper- 
iment: As an example, the calculated gap for NC's of 
30 A diameter is roughly the same as the experimental 
value 3.65 eV reported in Ref. [6j. It is interesting to note 
that, in the full size range, the HOMO and LUMO de- 
generacies we obtain are compatible with the T$ and Tq 
representations, respectively, in accord with the symme- 
try analysis of Ref. |28[ 

Another feature illustrated in Fig. Q] is a comparison of 
the electronic structure in two different crystallizations: 
Ga- and N-centered NC's. As can be seen in the figure, 
the difference between them is neglig ible for R > 15 A. 
Even in smaller NC's (see the insets) the small fluctua- 
tions do not seem to be significant. This behavior allows 
us to restrict our considerations to cation-centered NC's 
for an acceptor impurity and anion-centered NC's for a 
donor impurity even in the case of off-center impurities 
(see below) without any loss of generality. 

We next discuss the impurity states starting with the 
acceptor impurity. We calculate the HOMO in NC's with 
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FIG. 1: Highest valence band (HOMO) and lowest conduction 
band (LUMO) levels in Ga- and N-centered undoped NC's 
(open and closed symbols, respectively). In the insets the 
results for the smallest NC's are displayed. The dotted lines 
are drawn to guide the eye. 



a Mg atom replacing the Ga at the center and deduce the 
corresponding hole binding energy Efj as a function of the 
NC radius for different values of the parameter f/rj. The 
results are found to fit the simple expression (dropping 
the superscripts for simplicity) 



E b (R) = e b + A-R- 



(6) 



as shown in Fig. [3] for L/n = H cV. Note that this par- 
ticular value of Uq has been chosen for the Mg impu- 
rity in GaN, because the corresponding asymptotic limit 
E b (R — > oo ) = e b = 0.239 eV approximately reproduces 
the experimental bulk binding energy (0.236 eV). Com- 
plementary results concerning the QD acceptor states 
presented in Figs.|3l|4]and[8]correspond to the same value 
of Uq. The fitting parameters e b , A and /?, of course, de- 
pend on Uq. They may also depend on the QD size range. 
Generally speaking, the size ranges can be distinguished 
in terms of two different regimes of confinement: R > as 
(weak confinement) and R < as (strong confinement). In 
the strong-confinement regime the usual image of the im- 
purity center as a hydrogen-like entity is no longer valid. 
The asymptotic value -Eb(bulk) obviously corresponds to 
e b in the weak-confinement regime. Since the radii of all 
the investigated NC's are larger than the Mg acceptor 
Bohr radius, all of them belong to the weak-confinement 
regime in this case, and a single set of fitting parameters 
is sufficient as shown in Fig. 2. 

The degree of localization of the bound carriers is char- 
acterized through the analysis of wave functions and av- 
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FIG. 2: NC size-dependent acceptor binding energy for Uo = 
11 eV and the fitting curve E b = 0.239 + 23.149 R~ 1659 
(dashed line). 



erage radii. In Fig. [3] we show the shell- wise radial prob- 
ability distribution of the HOMO in an undoped QD of 
68 A radius in the upper panel. The ground-state of 
the hole in the same QD when doped with a Mg atom is 
shown in the lower panel. Both are band-edge wave func- 
tions, remarkably different from each other: The HOMO 
is spread over the whole NC with a maximum at r ~ -R/2, 
while the bound hole is essentially concentrated on the 
five nearest-neighbor atomic shells. This highly local- 
ized character of the acceptor hole in the ground state 
remains unchanged over the whole set of NC's, as illus- 
trated by the small size dependence of the average radius 
shown in Fig. 4, in agreement with experiment in the 
bulk. The magnetic resonance spectra in Mg-doped GaN 
layers are consistent with the same description of the 
hole, mainly distributed over the four nearest-neighbor 
atoms surrounding the impurity4^ The average radius of 
the hole orbit plotted as a function of the NC radius in 
Fig. 0] is found to fit the curve 



r a (R) 



Cl 



ci 



VEb(R) ^e b + C2-R- 



(7) 



where ci, ci and C3 are the fitting parameters. From 
([5]), this expression adequately reproduces the EM A bulk 
radius at large R: f a (R — > 00) = 6.3 A, which is indeed 
close to 3a%/2 = 6.7 A. 

The ground state of the donor impurity (O) is inves- 
tigated in a similar way. In order to obtain the value of 
Uq suitable for the O impurity in GaN we computed the 
binding energies for different values of Uq and fit their 
size dependence to the expression ([5]). Two distinct sets 
of fitting parameters are obtained for the two confine- 
ment regimes. In Fig. [5] we show the results for Uq = —2 
eV: The binding energy is strongly enhanced in the small- 
est NC's and tends to 0.027 eV in the bulk limit, which 
is close to the experimental value .E^ (bulk) = 0.033 eV. 
This value of Uq is then used in further analysis of the 
donor states: Figs. 6, 7 and 9. In Fig. [5] we show the 
probability distribution of the LUMO in an undoped 
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FIG. 3: Radial probability distribution of the HOMO in an 
undoped NC of 68 A radius (upper part), and of the ground- 
state acceptor hole in the same NC doped with a Mg atom 
(lower part), r is the radial distance from the QD center. 
Note that the two probability scales are different. Dotted 
lines guide the eye. 
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FIG. 4: QD size dependence of the average radius of the ac- 
ceptor hole in the ground state and the fitting curve f a = 
3.035/ V0. 235 + 31.685 R- 1 - 713 (dashed line). 
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FIG. 5: QD size-dependent donor binding energy for Uo = 
-2 eV and the fitting curves: E b = 0.078 + 22.386 R~ 1S37 in 
the strong-confinement regime (solid line) and Et = 0.027 + 
3.860 R~ ' in the weak-confinement (dashed line). 
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FIG. 6: Radial probability distribution of the LUMO in an 
undoped NC of 68 A radius (upper part) and that of the 
ground-state donor electron in the same NC when doped with 
an O atom (lower part). Note that the two probability scales 
are different. 



IV. RESULTS: OFF-CENTER IMPURITY 



QD of 68 A radius in the upper panel. That of the 
ground-state donor electron in the same QD when doped 
with an O atom is shown in the lower panel. They are 
the same band-edge states, but very different from each 
other, because of the carrier localization in the latter. 
From a comparison of Fig. [3] with Fig. [B] we see that 
the donor electron is much more spread out than the ac- 
ceptor hole, regardless of the crystal size. The average 
radius of the bound electron presented in Fig. [7] confirms 
this difference; it is remarkably larger than that of the 
bound hole in Fig. [4] for all the NC's. The fitting curve 
fd{R) in the weak-confinement regime tends to 40.7 A 
for R — > oo, in good agreement with the EMA radius in 
bulk, 3a|/2 = 41.6 A. 



We also compute off-center impurity states. As an 
example we focus on an intermediate-size NC of ra- 
dius 31.6 A, slightly larger than the donor Bohr radius 
(~ 28 A). The impurity position dependence of the ac- 
ceptor ground state is shown in Fig. 8, where the binding 
energy and the average radius are plotted against the dis- 
tance of the impurity atom from the NC center. Notice 
the slow evolution until the acceptor approaches the QD 
surface, where dramatic changes occur: The binding en- 
ergy drops and the average radius rises rather abruptly. 
This occurs when the Mg atom is placed within a surface 
layer of thickness ~ 4.5 A of the order of the acceptor 
Bohr radius. This behavior is related to the strongly 
localized nature of the bound hole state. 

On the other hand, in the case of donor impurity, the 
Bohr radius is comparable to the radius of the referred 
NC and the bound electron is spread over the whole NC. 
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FIG. 7: QD size-dependent average radius of the donor 
electron in the ground state and the fitting curves: f<j = 
6.697/V0.021 + 81.871 R- 1 - 912 in t he strong-confinement 
regime (solid line) and f d = 7.279/V0.032 + 81.495 i?" 1865 
in the weak-confinement (dashed line). 



As a result, the position dependence of the binding en- 
ergy and the average radius, as shown in Fig. [9j is quite 
smooth. The effects are important regardless of the O 
atom location. The closer the dopant gets to the crys- 
tal boundary, the smaller is the binding energy and the 
larger the average electron radius. This monotonously 
decreasing binding energy is comparable to that calcu- 
lated in Ref. l29l . From the experimental point of view, 
our results show that the acceptor-doped NC's are less 
sensitive to the impurity position than the donor-doped 
NC's. 
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FIG. 8: Binding energy and average radius of the hole ground 
state as a function of the acceptor position in a QD of radius 
31.6 A (closed and open symbols, respectively). The dashed 
and dotted lines are drawn to guide the eye. 
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FIG. 9: Binding energy and average radius of the electron 
ground state as a function of the donor position in a QD of 
radius 31.6 A (closed and open symbols, respectively). 



CONCLUDING REMARKS 



We have studied acceptor (Mg) and donor (O) states 
in GaN NC's doped with a single substitutional impurity 
within the sp 3 d 5 s* tight-binding model. The zinc-blende- 
structure crystallites are assumed spherical, ranging from 
1 to 13.5 nm in diameter. The computed binding en- 
ergy is highly enhanced with respect to the experimental 
bulk value when the dopant is placed at the center of the 
smallest QD's. For larger sizes it decreases following a 
scaling law that extrapolates to the bulk limit (236 meV 
for Mg and 33 meV for O). The degree of localization 
of the bound carriers is analyzed through their wave 
functions and average radii. The ground-state acceptor 
hole is mostly distributed over the nearest-neighbor anion 
shell. The donor electron, on the contrary, is much less 
localized. We also investigated off-center impurities in 
intermediate-size NC's (R ~ 32 A). The acceptor bind- 
ing energy is found to be barely dependent on the Mg 
impurity position unless it lies within a surface shell of 
~ 4.5 A thickness (~ the acceptor Bohr radius), where 
the ionization energy is drastically reduced. On the con- 
trary, the donor binding energy gradually decreases as 
the O dopant approaches the surface (the donor Bohr 
radius ~ 28 A is similar to the radius of the referred 
NC's). This difference arises from the larger spatial ex- 
tension of the bound electron as compared to the bound 
hole. Although we have focused on Mg and O impurities, 
our approach can be straightforwardly extended to other 
dopant species. 
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